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(54) Transform processing apparatus and method and medium for storing compressed digital signals. 

(5?) A forward transform processing method for 
performing a modified discrete cosine trans- 
form. A block of N time domain input samples is 
multiplied by a forward transform window to 
provide a block of N processed samples which 
is transformed into a block of N/2 spectral 
coefficients by preprocessing the block of N 
processed samples to provide a block of N/4 
complex numbers, transforming the block of 
N/4 complex numbers by the equivalent of a 
Fast Fourier Transform into a block of N/4 
transformed complex numbers that is then 
post-processed to provide an output block of 
N/2 spectral coefficients in the frequency do- 
main. The corresponding apparatus embodies 
the method just described. 

An inverse transform method for performing 
an inverse modified discrete cosine transform. 
A block of N/2 spectral coefficients is pre- 
processed to provide a block of N/4 frequency 
domain complex numbers, which is transfor- 
med into a block of N/4 time domain complex 
numbers by processing equivalent to a Fast 
Fourier Transform. The block of N/4 time do- 
main complex numbers is post processed to 
provide a block of N/2 processed samples in the 
time domain, which is re-arrayed and multiplied 
by an inverse transform window to provide an 
output block of N samples in the time domain. 
The apparatus embodies the method just des- 
cribed. 
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This invention relates to transform processing apparatus and methods and to a medium for storing a com- 
pressed digital signal. 

Orthogonal transform coding is known as a way of high efficiency coding sample data signals in the time 
domain, such as digital signals resulting from sampling analog audio signals. The input signal is divided in time 
into blocks of a predetermined number of samples. Transform coding orthogonally transforms each block of 
the input signal into a block of plural spectral coefficients in the frequency domain. Transform coding is typified 
by the discrete cosine transform (DCT). Discrete cosine transform coding suffers from distortion caused by 
discontinuous connecting points between adjacent blocks, and which manifests itself as noise. Block distortion 
has conventionally been reduced by overlapping one end of a block with the end of the adjoining block. A modi- 
fied discrete cosine transform (MDCT) overlaps a given block with the adjoining blocks on both sides by one 
half of the block length. The MDCT is well suited for high efficiency coding because samples in the overlapped 
portions are not transmitted. Thus, NI2 spectral coefficients are generated from a frame of N time domain sam- 
ples. 

Discussion of the MDCT and the inverse MDCT (IMDCT) may be found in Motizuki et al., Constraints in 
MDCT Filters in Which Plural Block Sizes Co-exist t SHINGAKU GIHOU, CAS90-10 and DSP90-14, pp. 55 to 
60, or in Sugiyama et al.. Adaptive Transform Coding with Adaptive Block Lengths Using MDCT (ATC-ABS), 
ABSTRACTS IN 1990 SPRING MEETING OF ELECTRONIC INFORMATION COMMUNICATION MEETING, 
A- 197. 

The modified discrete cosine transform and inverse modified discrete cosine transform will now be briefly 
described by referring to Figure 1. In Figure 1, a given block of sample data in the time domain, such as the 
J-tU block, has an overlap of one- half (50%) with each of the (J-1)-th block and the (J+1)-th block. If the number 
of samples in the fth block is N, where N is a natural number, there is an overlap of N/2 samples with the (J- 
1)-th block, and an overlap of NI2 samples with the )-th block. Each of the N samples in each block, such 
as the Mh block, is processed by a forward transform window or pre-filter Wh to produce a block 102 of N 
windowed samples in the time domain. 

The characteristics of the forward transform window or pre-filter Wh are so selected that power density 
of the transformed data is maximized in conformity with the statistical properties of input signals. The block 
102 of N windowed samples in the time domain in processed by a forward MDCT to produce a block 103 of 
independent spectral coefficients in the frequency domain. The number of spectral coefficients is one half the 
number of input samples in the block 102, or N/2. 

In the inverse transform, the block 103 of Ni2 spectral coefficients is processed with a linear IMDCT to 
produce a block 104 of N processed samples in the time domain. The block 104 of N processed samples is 
processed by an inverse transform window or synthesizing filter Wf to produce a block 105 of samples in the 
time domain. The block 105 of samples in the time domain is summed to the preceding and succeeding blocks 
of samples in the time domain to produce an output signal of samples in the time domain. 

A method of high speed processing using a Fast Fourier Transform (FFT) to provide an MDCT and IMDCT 
has been proposed by Iwatare et al., A Study in MDCT System and High Speed Processing Method, SHIN- 
GAKU-GIHOU, CAS 90-9, DSP 90-13, pp.49 to 54, and a similar method has also been proposed in Japanese 
Patent Early Publication 4-44099 (1992). 

The processing method proposed by Iwatare et al. will now be described. The MDCT is defined by the 
equation: 



Xq. MDCT input signal 
N: block length 

h: function of forward transform window 
yoi MDCT output signal 
Cq. constant, and 

n is an integer between 0 and N - 1, and k is an integer between 0 and N/2 - 1. The block numbers are omitted 
because MDCT processing is carried out on the time domain data in each block independently. The method 
for performing an MDCT is substantially independent of the value of so the following description assumes 
that C 0 = 1. 

In the calculation process, x 0 is multiplied by the forward transform window to find x 01 : 
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x 0 i(n) = x 0 (n)h(n) for O^n^N - 1 (2) 
x 02 is derived from x 01 in accordance with the following equation: 



4 4 



(3) 



Xq3 is calculated from in accordance with the following equation: 

*03(n) = Xo2(2n) - Xo 2 (A/ - 1 - 2n) for O^n^ - 1 (4) 
x 0 3 is multiplied by the coefficient shown in equation (5) to produce a block of complex numbers z 01 . 

(S) 



z 0 i is processed by a Fast Fourier Transform (FFT) having a length of N/2 to produce a block of complex 
numbers Zq 2 . 



/01 is a block of real numbers calculated from the block of complex numbers Zq 2 as follows: 



30 



y 01 thus calculated thus coincides with y 0 as defined in equation (1), as demonstrated by the above-cited 
35 treatise A Study in MDCT System and High Speed Processing Method. 
Conversely, the IMDCT is defined by the following equation (8): 



where 

yv IMDCT input signal 
45 N: block length 

t inverse transform window function 
x t : IMDCT output signal 
C A : constant, and 

where n is an integer between 0 and N - 1, and k is an integer between 0 and N/2 - 1 . The block numbers are 
so omitted because IMDCT processing is carried out on the frequency domain data in each block independently. 
The present IMDCT calculating method is substantially independent of the value of C 1( so the following de- 
scription assumes that = 1. 

In the known calculation method, y 1 is re-arrayed in accordance with the following equation to produce 

yn- 

55 
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(9) 



y u is multiplied by the coefficient shown in equation (10) to produce a block of complex numbers z u . 
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*n<*)=*i<*)«p(-^) for 0***^-1 



(10) 



15 z u is processed by an inverse FFT of a length of N/2 to produce a block of complex numbers z 12 . 



20 



(11) 



The block of complex numbers z 12 is multiplied by the coefficient shown in equation (12), and the real part 
is extracted to produce the block of real numbers x n . 



25 



tll (n^Re\2z l2 (n)tx P \-j^^ for 0 



(12) 



30 x,, is re-arrayed with a sign change and multiplied by the inverse transform window in accordance with 

the equation: 
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This coincides with x, as defined in equation (8), as demonstrated by the above-cited treatise A Study in 
MDCT System and High Speed Processing Method. 

The processing to perform an MDCT to determine the spectral coefficients for block of N samples and to 
perform a complementary IMDCT involves performing a complex number FFT having a length of N/2. Suitable 
algorithms for this are described in E. Oran Brigham, THE FAST FOURIER TRANSFORM. In chapter 11, en- 
titled Theory of FFT Algorithm with a Base of 2, the algorithm described on pages 196 to 199 in the Japanese 
translation by Miyagawa et al. has the following requirements: 
/V74 log(A//2) complex number multiplications, 
N/2 log(A//2) complex number additions, and 
a work area for storing N/2 complex numbers. 

Thus, the method described above is not fully suitable for high-speed processing, especially when a small 
processing system is used. 

In accordance with a first aspect of the present invention, a forward transform processing apparatus is 
provided for carrying out a modified discrete cosine transform on a block of N samples of a digital input signal 
in the time domain to provide a block of N/2 spectral coefficients in the frequency domain. The apparatus com- 
prises a circuit that multiplies the block of N samples by a forward transform window to provide a block of N 
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forward transform window processed samples, and a transform circuit that transforms the block of N forward 
transform window processed samples to provide a block of N/2 spectral coefficients. The transform circuit in- 
cludes a preprocessing circuit that preprocesses the block of N forward transform window processed samples 
to provide a block of N/4 complex numbers. The transform circuit also includes an FFT circuit that performs 

5 processing equivalent to a Fast Fourier Transform on the block of N/4 complex numbers to provide a trans- 
formed block of N/4 complex numbers. Finally, the transform circuit includes a postprocessing circuit that post- 
processes the transformed block of complex numbers to provide the block of N/2 spectral coefficients. 

In accordance with a second aspect of the present invention, an inverse transform processing apparatus 
is provided for carrying out a inverse modified discrete transform on a block of N/2 spectral coefficients in the 

10 frequency domain to provide a block of N samples in the time domain as a digital output signal. The apparatus 
comprises an inverse transform circuit that inversely transforms the block of N/2 spectral coefficients to provide 
a block of N/2 processed samples in the time domain, and a circuit that re-arrays and multiplies the block of 
N/2 processed samples in the time domain by an inverse transform window to provide the block of N samples 
in the time domain. The inverse transform circuit includes a pre-processing circuit that pre-processes the s pee- 
rs tral coefficients in the block of N/2 spectral coefficients to provide a block of N/4 frequency domain complex 
numbers, and an FFT circuit that performs processing equivalent to a Fast Fourier Transform on the block of 
N/4 frequency domain complex numbers to provide a block of N/4 time domain complex numbers. Finally, the 
inverse transform circuit includes a post-processing circuit that post- processes the block of N/4 time domain 
complex numbers to provide the block of N/2 processed samples in the time domain. 

20 In accordance with a third aspect of the present invention, a method is provided for performing a modified 

discrete cosine transform on a block of N samples of a digital input signal in the time domain to provide a block 
of N/2 spectral coefficients in the frequency domain. According to the method, the block of N samples is mul- 
tiplied by a forward transform window to provide a block of N forward transform window processed samples, 
and the block of Nforward transform window processed samples is transformed to provide a block of N/2 spec- 

25 tral coefficients. The step of transforming the block of Nforward transform window processed samples includes 
the steps of preprocessing the block of N forward transform window processed samples to provide a block of 
N/4 complex numbers; performing processing equivalent to a Fast Fourier Transform on the block of N/4 com- 
plex numbers to provide a transformed block of N/4 complex numbers; and post-processing the transformed 
block of complex numbers to provide the block of N/2 spectral coefficients. 

30 In accordance with a fourth aspect of the present invention, a method is provided for performing an inverse 

modified discrete transform on a block of N/2 spectral coefficients in the frequency domain to provide a block 
of N samples in the time domain as a digital output signal. According to the method, the block of N/2 spectral 
coefficients is inversely transformed to provide a block of N/2 processed samples in the time domain, and the 
block of N/2 processed samples in the time domain is re-arrayed and multiplied by an inverse transform window 

35 to provide the block of N samples in the time domain. The step of inversely transforming the block of N/2 spec- 
tral coefficients includes the steps of pre-processing the spectral coefficients in the block of N/2 spectral coef- 
ficients to provide a block of N/4 frequency domain complex numbers, performing processing equivalent to a 
fast Fourier Transform on the block of N/4 frequency domain complex numbers to provide a block of N/4 time 
domain complex numbers, and post-processing the block of N/4 time domain complex numbers to provide the 

40 block of N/2 processed samples in the time domain. 

In accordance with a fifth aspect of the present invention, a medium is provided for storing a compressed 
digital signal comprising a block of N/2 spectral coefficients. The block of N/2 spectral coefficients is derived 
from a block of N samples of a digital input signal in the time domain by a modified discrete cosine transform 
performed according to a method that includes the steps of multiplying the block of N samples by a forward 

45 transform window to provide a block of N forward transform window processed samples, and transforming the 
block of N forward transform window processed samples to provide a block of N/2 spectral coefficients. The 
step of transforming the block of N forward transform window processed samples includes the steps of pre- 
processing the block of N forward transform window processed samples to provide a block of N/4 complex 
numbers; performing processing equivalent to a Fast Fourier Transform on the block of N/4 complex numbers 

50 to provide a transformed block of N/4 complex numbers; and post- processing the transformed block of complex 
numbers to provide the block of N/2 spectral coefficients. 

In the method and apparatus according to the invention, a linear forward transform for providing an MDCT 
of N samples of a digital input signal may be realized by applying FFT processing to complex signals each hav- 
ing a length of N/4. This significantly reduces the number of processing operations and the work area required. 

55 Additionally, in the apparatus and method according to the invention, a linear inverse transform for pro- 

viding an IMDCT of N/2 independent spectral coefficients may be performed by applying FFT processing to 
complex signals each having a length of N/4. This also significantly reduces the number of processing oper- 
ations and work area required. 

5 
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The invention will be described hereinafter, by way of example only, with reference to the accompanying 
drawings, in which: 

Figure 1 is a diagrammatic view showing the processing steps for performing a modified discrete cosine 
transform (MDCT) and an inverse modified discrete cosine transform (IMDCT). 
5 Figure 2 is a flow chart schematically showing the processing steps of a method of performing an MDCT 

according to the present invention. 

Figure 3 is a functional block diagram of an apparatus according to the invention for implementing the 
method of performing an MDCT shown in Figure 2. 

Figure 4 is a flow chart schematically showing the processing steps of a method of performing an IMDCT 
io according to the present invention. 

Figure 5 is a functional block diagram of an apparatus according to the invention for implementing the 
method of performing an IMDCT shown in Figure 4. 

Figure 6 is a block diagram showing an example of a circuit arrangement of a high efficiency encoding 
apparatus to which an apparatus and method for performing an MDCT according to the invention is applied. 
15 Figure 7 is a block diagram showing an example of a circuit arrangement of a high efficiency decoding 

apparatus to which an apparatus and method for performing an IMDCT according to the invention is applied. 

The calculation principles employed in the present invention will be described before practical embodi- 
ments of the apparatus and the method according to the present invention are described. 

To calculate an MDCT using an apparatus and a method according to the present invention, a block of 
20 N/4 complex numbers Z 0 i, given by the following equation (14), 

Z ox {n) -(x 03 (2n) + jx 02 (2ml) ) exp (-j-JL*£_) for 0sm|-1 11 ^ 

25 

is derived from x 03 of equation (4), and a block of N/2 complex numbers Zq 2 is calculated in accordance with 
the equation: 

n/a - 1 

30 Z„(k). £ Z 0l (*) exp(-j-i™£L) for OsJcs|-l (15) 

From the block of complex numbers Zo 2 of equation (15), a block of complex numbers Z03 and a block of 
complex numbers Z04, each block having a length N/2, are calculated in the following manner 

35 / 

W-^OO^j-l-*)) M OikiZ-l (16) 



^.(l^^l^-^fl-*))]^ 0 «nfl (17) 

Calculating a block of N/2 real numbers y 01 from Z03 and Z04, we obtain: 

Vt)^^^^^^^*^*))) for 0<**|-1 (18) 

This coincides with y 0 defined by equation (1). This may be seen from the following: by substituting equa- 
tion (14) in equation (15), we get 
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• 1' (%c»> •*«*•*')) «p(-/ a ' ,Q * ) ;" 2 " ) (19) 



and 

10 



- £ fa^)-^^^))^^ 21 ^ 1 ^ ) (20) 



r5 



20 



25 



30 



Substituting equations (19) and (20) in equations (16) and (17), we obtain for Zq 3 and 2^: 
Z„<0 = £ % (2»)e^(-;^^p2) (2D 



Subtituting the values of Zq Z and Z04 into equation (18), we obtain for YTJ1: 

.*(1'>„(»)«p(-/2=«!11>| 



35 / 
/ 



40 



45 



From equations (1) through (4) we obtain: 
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2*Qk+ti 



•14-f) 

• g' (*»c2i.)-^w-i-a.))«»(250*^pii)J 

■g^H ^y^ ) (24) 



Equations (23) and (24) clarify that Y 0 i coincides with y 0 and show that an MDCT may be performed using 
a method according to the present invention. 

40 Since Z 02 has a period of N/4, the calculating process may be performed for values of k from 0 to (N/4) • 

1. An FFT having a length of N/4 may be used for this. Additionally, a shown in equation (21), since Z 03 has a 
period of N/4, the process for calculating equation (16) may simply be carried out for values of k from 0 to (N/4) 
- 1. Similarly, as shown in equation (22), since has a period in which only the sign is changed at an interval 
of N/4, the process for calculating equation (1 7) may simply be carried out for values of k from 0 to (N/4) - 1 . 

45 These simplifications enable the calculating process to be carried out more efficiently. 

Further, when performing an I MDCT, the block of N/4 complex numbers Z n given by the following equa- 
tion: 



50 



55 



"(y„(2*) *A u (2*M))e«p(-y^) for 0<*i£-l (25) 

is determined form y n of equation (9), and the block of N/2 complex numbers Z 12 is calculated in accordance 
with equation (26): 
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(20 



From the block Z 12 . a block of A//2 complex numbers Z 13 and a block of Nil complex numbers Z H are cal- 
culated using equations (27) and (28) as follows: 



10 



15 



20 



W 9 {( z u<«) + ^7-i-")) o*«£-i 



(27) 



1 (28) 



The block of N/2 real numbers is calculated from Z 13 and Z 14 in accordance with: 

X„(/i) = fc (wp(-/ a "^ l} ) (Z|»W^ M W)) A* (29) 



The block of real numbers X u thus calculated is processed by sign exchange, re-arraying and inverse 
25 transform window application to calculate the block of N real numbers X n to obtain: 
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for 0S«£-1 



(30) 



Equation (30) shows that Y 12 coincides with x 1 defined by equation (8) by the following analysis. Substi- 
tuting equation (25) in equation (26), we obtain the following equations: 



(Jl) 



50 



55 



i--) • E (yn^-^^i))^ ^ 2 ^-^^) 2 * 



(32) 
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Substituting equations (33) and (34) into equation (29), and using equation (9), we obtain: 
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/ 



50 Applying the relation defined by equation (30) to equation (35), we obtain equations (36), (37), and (38), 

depending on the value of n: 
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•f(n) £ y^k) cos 



Jf-O 



nl2k*l) (2 



2N 



to 



15 



jr/2-i ( n(2iol) (2(/3--^) + l) + * (2A>1)2N 
f<n> E yi<«cos 1 ±1 



Jc-0 



N/2-1 



2N 



- f (n) g y x (W cos - 



n (2Jc+l) 



(an.fi)] 



2tf 



for ^zn*N-l 
(38) 



20 



Comparing equations (36), (37), and (38) with equation (8) shows that an IMDCT may be performed using 
the method of the present invention. 

Since Z 12 has a period of N/4, as shown in equation (26), the above calculating process may simply be 
performed for values of n from 0 to N/4 - 1. An FFT having a length of N/4 may be used for this. Additionally, 
as shown in equation (33), since Z 13 has a period of N/4, the process for calculating equation (27) may simply 
be carried out for values of n from 0 to N/4 - 1 . Similarly, as shown in equation (34), since Z u has a period in 
which only the sign is changed at an interval of N/4, the process for calculating equation (28) may simply be 
25 carried out for values of n from 0 to N/4 - 1 . These simplifications enable the calculating process to be carried 
out more efficiently. 

It can be seen that, with an apparatus and the method according to the present invention, the MDCT de- 
fined by the equation (1) and the IMDCT defined by the equation (3) may be calculated using an FFT having 
a length of N/4. Such an FFT may be executed by 
N/8 log (N/4) complex number multiplications, 
N/4 log (N/4) complex number additions, and 
a work area for storing N/4 complex numbers. 
This means that the number of the processing operations and the work area may be significantly reduced 
relative to the conventional method. 

In performing the MDCT, the operations for calculating equations (16) to (18) may be combined into equa- 
tion (39): 



30 



35 



40 



45 



Y ox lk) 



10* (2*+l) 
AN 



)) a (MI- l -1 



(39) 



50 



to find y 01 directly from Z 02 . 
Additionally, if equation (40): 



55 



ziz(k) - g z 01 (n) exp(jl2JSH) for 0 



(40) 



is calculated for an inverse FFT instead of calculating equation (15) using a forward FFT, Y 0 i may be found 
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from Zw, because Z$i can be obtained by re-arraying Z Q2 . 

A method according to the present invention provides the processing for performing an MDCT with an 

input of N samples using processing equivalent to an FFT having a length of A//4, together with pre-processing 

and post-processing. Performing the MDCT using equations (39) and (40) is therefore encompassed by the 
5 present invention. The same may be said of the processing operations for performing the IMDCT. Thus, the 

method of executing the processing operations for calculating equations (27) through (30), or of using an I FFT 

in place of the FFT of equation (26) is encompassed by the present invention. 

Referring to the drawings, a preferred embodiment of the method according to present invention will next 

be described. Figure 2 is a flow chart schematically illustrating the sequence of signal transform operations 
w in a method for performing an MDCT according to the present invention. 

In the first step S01 of Figure 2, sample data in the time domain, such as PCM audio data, is divided into 

blocks. Each block consists of a predetermined number of samples. In the example shown in Figure 2, each 

block consists of N samples. The blocks are arranged such that there is a 50% overlap between adjacent 

blocks, i.e., adjacent blocks are overlapped by N/2 samples, as shown in Figure 1. The block of N samples x 0 
15 is multiplied by the forward transform window Wh, as shown in Figure 1, to provide a block of N windowed 

samples Xqv 

In the next step, S02, the block of N windowed samples x 01 is converted by the transform defined by equa- 
tion (4) into a block of N/2 real numbers x 03 . 

In the next step S03, the even-numbered elements and the odd-numbered elements of x 03 are allocated 
20 to real and imaginary parts, respectively, of a block of N/4 complex numbers. The resulting complex numbers 
are then multiplied by the complex coefficient shown in equation (14) to convert them into the block of N/4 
complex numbers 2^. 

In the next step S04, Z01 is processed by an FFT to produce the block N/4 complex numbers in the 
frequency domain. 

25 In the next step S05, the block of complex numbers 2o 2 is processed in accordance with equations (16) 

through (18) to provide the block of N/2 spectral coefficients Vqi, which are the output of the MDCT. 

Figure 3 shows a typical hardware arrangement for carrying out the method of performing an MDCT ac- 
cording to the invention. A block of sample data in the time domain stored in a time domain sample buffer is 
converted into spectral coefficients in the frequency domain by the MDCT processing circuit. 
30 The circuit shown in figure 3 can be realized using discrete circuitry, or multiple microprocessors. Prefer- 

ably, however, the majority of the processing operations are performed using an appropriately programmed 
digital signal processor circuit and auxiliary memories. 

In Figure 3, samples in the time domain, such as a PCM audio signal, are supplied to the input terminal 
11, whence they pass into the time domain sample buffer 12. In the time domain sample buffer, the samples 
35 are divided into blocks of a predetermined number of samples, for example, into blocks of N samples. The 
blocks of N samples are then supplied to the MDCT processing circuit 20. 

The MDCT processing circuit 20 includes an equation (2) calculating circuit 21, an equation (3) calculating 
circuit 22, an equation (4) calculating circuit 23, an equation (14) calculating circuit 24, and an equation (15) 
calculating circuit 25 for calculating the equations (2), (3), (4), (14), and (15), respectively. The calculated re- 
40 suits from the equation (16) calculating circuit 25 are supplied to both the equation (16) calculating circuit 26 
and to the equation (17) calculating circuit 27. Outputs Zq 3 from the equation (16) calculating circuit 26 and 
Z04 from the equation (17) calculating circuit 27 are supplied to the equation (18) calculating circuit 28, which 
calculates Y 0h and provides V 0 i to the output terminal 13. 

Figure 4 is a flow chart schematically showing the sequence of signal processing for carrying out a method 
45 of performing an IMDCT according to the present invention. 

In the first step S11 in Figure 4, a block of N/2 spectral coefficients y 1 are processed by sign changing 
and re-arraying, as shown in equation (9), and are converted into the block ofA//2 real numbers y u . 

In the next step S12, the even-numbered elements and the odd-numbered elements of y n are allocated 
to real and imaginary parts, respectively, of a block of N/4 complex numbers. The block of N/4 complex numbers 
so is then multiplied by the complex coefficient shown in equation (25) to provide the block of N/4 complex num- 
bers Zn. 

In the next step S13, Z n is processed by an FFT to produce the block of N/4 complex numbers Zi 2 in the 
frequency domain. 

In the next step S14, Z 12 is processed in accordance with equations (27) through (29) to produce the block 
55 of N/2 real numbers X n . 

In the next step S1 5, the block of N/2 real numbers X, t is processed by sign change, rear raying, and inverse 
transform window multiplication in accordance with equation (30) to produce the block of N real numbers X 12 , 
from which samples of an output signal in the time domain can be derived. 

14 
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Figure 5 shows a topical hardware arrangement for carrying out a method of performing an IMDCT trans- 
form according to the invention, and for constructing a digital output signal in the time domain. 

The circuit shown in Figure 5 can be realized using discrete circuitry, or multiple microprocessors. Prefer- 
ably, however, the majority of the processing operations are performed using an appropriately programmed 
5 digital signal processor circuit and auxiliary memories. The digital signal processor used for the circuit shown 
in Figure 3 can be reprogrammed to provide the processing operations of the circuit shown in Figure 5. 

In the circuit shown in Figure 5, the block of NI2 spectral coefficients /i, resulting from an MDCT, is supplied 
to the input terminal 16. Normally, the block of spectral coefficients will have been derived from a com- 
pressed digital signal reproduced from a storage medium, such as a magneto-optical disk, or from a transmis- 
w sion medium. 

From the input terminal, the block of spectral coefficients yi is supplied to the equation (9) calculating cir- 
cuit 31 of the IMDCT calculating circuit 30. The output of the equation (9) calculating circuit 31 is supplied via 
the equation (25) calculating circuit 32 to the equation (26) calculating circuit 33. The output Z 12 from the equa- 
tion (26) calculating circuit 33 is supplied to both the equation (27) calculating circuit 34 and the equation (27) 

75 calculating circuit 35. The outputs from the equation (27) calculating circuit 34 and the equation (28) calculating 
circuit 35 are supplied to the equation (29) calculating circuit 36 which provides an output X u to the equation 
(30) calculating circuit 37. The output, X 12 , from the equation (30) calculating circuit 37 is provided as the output 
of the IMDCT calculating circuit 30. 

The output of the IMDCT calculating circuit 30, the block of N real numbers X 12 is supplied to the overlap 

20 and add circuit 18 via the time domain sample buffer 17. Adjacent blocks of samples in the time domain data 
are overlapped by N/2 samples, and are added together to provide an output signal consisting of samples in 
the time domain. 

Figure 6 shows a practical example of a high efficiency coding apparatus in which the method and appa- 
ratus for performing an MDCT is applied. The high efficiency coding apparatus shown in Figure 6 uses the 

25 techniques of sub-band coding, adaptive transform coding, and adaptive bit allocation. 

A digital input signal, such as a digital PCM audio signal, is divided in frequency into a frequency range 
signal in each of plural frequency ranges. The bandwidths of the frequency ranges broaden in the direction of 
increasing frequency. An MDCT is performed on the frequency range signal in each of the frequency ranges 
independently. The resulting spectral coefficients in the frequency domain are grouped into critical bands, and 

30 the spectral coefficients in each critical band are coded by adaptive bit allocation. 

In Figure 6, an audio pCM input signal in the frequency range of, e.g., 0 Hz to 20 kHz are supplied to the 
input terminal 51 . The band dividing filter 52, which is preferably a Quadrature Mirror Filter (QMF), divides the 
input signal in frequency into a 0 Hz to 10 kHz frequency range, and a 10 to 20 kHz frequency range. A second 
QMF 53 divides the 0 Hz to 10 kHz frequency range in frequency into a 0 Hz to 5 kHz frequency range and a 

35 5 to 10 kHz frequency range. The 10 to 20 kHz frequency range signal is transmitted from the QMF 52 to the 
modified discrete cosine transform circuit (MDCT circuit) 54. The 5 to 10 kHz frequency range signal from the 
QMF 53 is transmitted to the MDCT circuit 55, and the 0 Hz to 5 kHz frequency range signal from the QMF 53 
is transmitted to the MDCT circuit 56. The MDCT circuits 54, 55, and 56 embody examples of a method and 
apparatus according to the invention. 

40 The spectral coefficients in the frequency domain, resulting from orthogonally transforming the three fre- 

quency range signals using the MDCT circuits 54, 55, and 56, are grouped into critical bands and supplied to 
the adaptive bit allocation circuit 57. Grouping the spectral coefficients into critical bands enables the masking 
characteristics of the human auditory sense to be taken into account 

Bits are allocated to the spectral coefficients in a way that takes advantage of a psychoacoustic property 

45 of the human auditory sense called "masking." Masking is a psychoacoustic phenomenon in which a sound is 
rendered inaudible, or "masked," by other sounds occuring simultaneously with, or slightly earlier than, or later 
than, the sound. Masking effects may be classed into time domain masking effects, that is, masking by sounds 
occurring earlier or later than the masked sound, and concurrent masking effects, which is masking is by si- 
multaneously-occurring sounds having a frequency different from the frequency of the masked sound. 

so Masking enables a sound to render inaudible any noise within its time or frequency masking range. This 

means that, in the presence of a signal that, when reproduced, produces a sound, a digital coding system may 
produce a quantizing noise level that is high compared with the noise level that is allowable in the absence 
of the signal, provided that the quantizing noise lies within the masking range of the sound produced by the 
signal. Since relatively high levels of quantizing noise are allowable if masked by the sound resulting from the 

55 signal, the number of bits required to quantize the signal representing the sound, or parts of the signal, may 
be significantly reduced. 

A critical band is a frequency band that takes advantage of the masking characteristics of the human au- 
ditory sense. A critical band is the band of noise that can be masked by a pure sound that has the same intensity 
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as the noise and has a frequency in the vicinity of the frequency of the noise. The width of the critical band 
increases with increasing frequency of the pure sound. The entire audio frequency range of 0 Hz to 20 kHz 
can be devided into, for example, 25 critical bands. 

In the adaptive bit allocation coding circuit 57, the spectral coefficients are first normalized by applying 

5 block floating to the spectral coefficients in each critical band. A block floating coefficient e.g., the maximum 
of the absolute values of the spectral coefficients in the critical band, is determined for each critical band. The 
spectral coefficients in each critical band are then normalized using the block floating coefficient. The nor- 
malized spectral coefficients in each critical band are then re-quantized using a number of bits that is large 
enough to provide a quantizing noise level that is low enough to be masked by the signals in the critical band 

10 and neighbouring critical bands. The output signal provided to the output terminal 58 comprises the adapt ively 
quantized spectral coefficients, together with the block floating coefficient for each critical band, and word 
length data indicating the number of bits used for coding the spectral coefficients in each critical band. The 
output signal can then be transmitted or recorded on a suitable recording medium for later reproduction and 
decoding. Examples of suitable recording media are magneto-optical disks, optical disks, magnetic tapes, and 

is semiconductor memories. 

Figure 7 shows a decoding apparatus complementary to the high efficiency coding apparatus shown in 
Figure 6. The decoding apparatus includes IMDCT circuits to which the IMDCT apparatus and method accord- 
ing to the invention is applied. Compressed digital data, which is normally reproduced from a storage medium, 
such as a magneto-optical disk, an optical disk, a magnetic tape, or a semiconductor memory, or a transmission 

20 medium, is fed into the terminal 61 . The compressed digital data includes a block floating coefficient for each 
critical band, word length data including the number of bits used for coding each critical band, and adaptively 
quantized spectral coefficients. Blocks of spectral coefficients are constructed from this data by a spectral 
coefficient decoding circuit 62. 

The spectral coefficients in each block are grouped into the spectral coefficients for the 1 0 to 20 kHz fre- 

2$ quency range, the spectral coefficients for the 5 to 10 kHz frequency range, and the spectral coefficients for 
the 0 Hz to 5 kHz frequency range for feeding into the IMDCT circuits 63, 64 and 65, respectively. The IMDCT 
circuits 63, 64, and 65, which employ an IMDCT processing method and apparatus according to the invention, 
transform the groups of spectral coefficients in the frequency domain into three frequency range signals in 
the time domain. 

30 Of the frequency range signals obtained as just described, the 0 Hz to 5 kHz frequency range signal and 

the 5 to 10 kHz frequency range signal waveform are first combined by the band combining circuit 66 into a 
0 Hz to 10 kHz frequency range signal. The 0 Hz to 10 kHz frequency range signal is combined with the 10 to 
20 kHz frequency range signal in the band combining circuit 67. The resulting full frequency range signal is 
provided supplied to the output terminal 68. 

35 it is to be noted that the present invention is not limited to the above described embodiments. For example, 

the present invention is not limited to use in the high efficiency coding and decoding apparatus shown in Fig- 
ures 6 and 7, but may be applied to a variety of transform coding apparatus or decoding devices. 

It is seen from above that the present invention can provide a high efficiency coding apparatus for audio 
signals in which a block of N time domain samples data overlapped 50% with a preceding and a subsequent 

40 block may be efficiently orthogonally transformed using a forward modified DCT (MDCT) which uses a Fast 
Fourier Transform (FFT) with a length of A//4. This significantly reduces the number of processing operations 
and the work area required. 

The present invention can also provide a decoding apparatus and method in which N/2 independent spec- 
tral coefficients are processed by an inverse modified DCT (MDCT) which uses a FFT with a length ofA//4. 

45 This also significantly reduces the number of processing operations and the work area required. 

Thus, an embodiment of the invention provides a signal processing apparatus and method wherein the 
processing operations for performing an MDCT and an IMDCT are carried out using an FFT with a shorter 
length. This reduces the number of processing operations, such as the number of multiplications and the num- 
ber of additions, reduces the work area required, and achieves a higher operating speed. 

50 

Claims 

1. A forward transform processing apparatus for carrying out a modified discrete cosine transform in a block 
55 of N samples of a digital input signal in the time domain to provide a block of NI2 spectral coefficients in 

the frequency domain, the apparatus comprising: 

(1) a means for multiplying the block of N samples by a forward transform window to provide a block 
of N forward transform window processed samples; and 
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(2) a transform means for transforming the block of N forward transform window processed samples 
to provide a block of N/2 spectral coefficients, the transform means including: 

(a) a preprocessing means for preprocessing the block of N forward transform window processed 
samples to provide a block of N/4 complex numbers, 

(b) an FFT means for performing processing equivalent to a Fast Fourier Transform on the block 
of N/4 complex numbers to provide a transformed block of N/4 complex numbers; and 

(c) a postprocessing means for post- processing the transformed block of complex numbers to pro- 
vide the block of N/2 spectral coefficients. 

The forward transform processing apparatus of claim 1 ? wherein, 

(1 ) the block of N forward transform window processed samples includes an n-th sample, a (n + 3A//4)- 
th sample, the (n + 3/v74)-tn samples having a sign, and a (n - N/4)-th sample; and 

(2) the preprocessing means includes: 

(a) a first deriving means for deriving, from the block of N forward transform window processed sam- 
ples, a block of first intermediate samples, the block of first intermediate samples including an n- 
th a (N - 1 - 2n)-th, and a 2n-th first intermediate sample, the first deriving means deriving the block 
of first intermediate samples by: 

(i) reversing the sign of the (n + 3/v74)-th sample to provide a sign-reversed (n + 3A//4)-th sample, 
and providing the sign-reversed (n + 3A//4)-th samples as the n-th first intermediate sample, for 
each integral value of n from 0 to N/4 - 1 , and 

(ii) providing the (n + A//4)-th samples as the nth first intermediate sample, for each integral value 
of n from N/4 to N - 1 , 

(b) a second deriving means for deriving, from the block of first intermediate samples, a block of 
second intermediate samples, the block of second intermediate samples including an n-th, a 2n- 
th, and a (2n + 1)-th second intermediate sample, the second deriving means deriving the block of 
second intermediate samples by subtracting, for each integral value of n from 0 to N/2 - 1 ,' the (N - 

1 - 2n)-th first intermediate sample from the 2n-th first intermediate sample to provide the n-th 
second intermediate sample; 

(c) a third deriving means for deriving, from the block of second intermediate samples, the block of 
N/4 complex numbers, the complex numbers in the block of N/4 complex numbers each having a 
real part and an imaginary part, and including an n-th complex number, the third deriving means 
deriving the block of N/4 complex numbers by: 

(i) selecting the 2n-th second intermediate samples as a real part of an n-th intermediate com- 
plex number, 

(ii) selecting the (2n + 1)-th second intermediate sample as an imaginary part of the n-th inter- 
mediate complex number, and 

(iii) multiplying the n-th intermediate complex number by an exponential of (-4njnfN) to provide 
the n-th complex number in the block of N/4 complex numbers, for each integral values of n from 

. OtoW/4-1. 

The forward transform processing apparatus of claim 1 , wherein, 

(1) the block of N/2 spectral coefficients includes a /c-th spectral coefficient, where k is an integer be- 
tween 0 and N/2 - 1, and 

(2) the post- processing means includes: 

(a) a means for expanding the transformed block of N/4 complex numbers based on periodically to 
provide a block of 3A//4 complex numbers, the block of 3A//4 complex numbers including a /c-th com- 
plex number and a (N/4 - 1 - /c)-th complex number, 

(b) a means for deriving a (N/4 - 1 - /c)-th complex conjugate from the (N/4 -1- /c)-th complex number, 

(c) a means for 

(i) summing the (N/4 - 1 - /c)-th complex conjugate and the /c-th complex number to provide a 
sum, and 

(ii) dividing the sum by two to provide a /c-th third intermediate complex number, 

(d) a means for 

(i) subtracting the (N/4 - 1 - /c)-th complex conjugate from the /c-th complex number to provide 
a difference, 

(it) dividing the difference by 2j to provide a quotient, where j is the square root of minus one, 
and 

(iii) multiplying the quotient by the exponential of (-2nj(2k + 1)/W) to provide a /c-th fourth inter- 
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mediate complex number, and 
(e) a means for 

(i) summing the k-th third intermediate complex number and the k-th fourth intermediate com- 
plex number to provide a k-th fifth intermediate complex number, and 

(ii) multiplying the k-th fifth intermediate complex number by an exponential of {-nj{2k + 1)/(2A/)) 
to provide a Ac-th sixth intermediate complex number having a real part and an imaginary part, 
and 

(iii) providing the real part of the k-th sixth intermediate complex number sum as the Ar-th spec- 
tral coefficient. 

The forward transform processing apparatus of claim 1, wherein 

(1) the block of N/2 spectral coefficient includes a k-th spectral coefficient, where k is an integer be- 
tween 0 and N/2 - 1, and 

(2) the post-processing means includes: 

(a) a means for expanding the transformed block of N/4 complex numbers based on periodicity to 
provide a block of 3N/4 complex numbers, the complex numbers in the block of 3N/4 complex num- 
bers having a real part and an imaginary part, and including a k-th complex number and a (N/4 - 1 
- k)-th complex number 

(b) a means for determining 

(i) a sum, the sum being one half of a sum of the cosine of (2n(2k + 1)/4N) and sine of (10*(2* 
+ 1)/4N), 

(ii) a first difference, the first difference being one half of a difference between the cosine of 
(2iz(2k + 1)/4N) and the sine of (10k(2* + 1)/4N), and 

(iii) a second difference, the second difference being one half of a difference between and sine 
of (10n(2Ac+ 1)/4N) and the cosine of (2*(2* + 1)/4N); 

(c) a means for multiplying: 

(i) the real part of the k-th complex number by the first difference to provide a first intermediate 
number, 

(ii) the real part of the (N/4 - 1 -*)-th complex number by the sum, to provide a second inter- 
mediate number, 

(iii) the imaginary part of the k-th complex number by the sum to provide a third intermediate 
number, and 

(iv) the imaginary part of the (N/4 - 1 - k)-th complex number by the second difference to provide 
a fourth intermediate number; and 

(d) a means for summing the first, the second, the third, and the fourth intermediate numbers 
to provide the Ac-th spectral coefficient. 

An inverse transform processing apparatus for carrying out a inverse modified discrete cosine transform 
on a block of N/2 spectral coefficients in the frequency domain to provide a block of N samples in the 
time domain as a digital output signal, the apparatus comprising: 

(1) an inverse transform means for inversely transforming the block of N/2 spectral coefficients to pro- 
vide a block of N/2 processed samples in the time domain, the inverse transform means including: 

(a) a pre-processing means for pre-processing the spectral coefficients in the block of N/2 spectral 
coefficients to provide a block of N/4 frequency domain complex numbers, 

(b) an FFT means for performing processing equivalent to a Fast Fourier Transform on the block 
of N/4 frequency domain complex numbers to provide a block of N/4 time domain complex numbers, 
and 

(c) a post- processing means for post-processing the block of N/4 time domain complex numbers 
to provide the block of N/2 processed samples in the time domain; and 

(2) a means for re-arraying and multiplying the block of N/2 processed samples in the time domain by 
an inverse transform window to provide the block of N samples in the time domain. 

The inverse transform processing apparatus of claim 5, wherein, 

(1) the block of N/2 spectral coefficients includes a 2/r-th, and a (N - 1 - 2*)-th spectral coefficient, 
each having a sign, 

(2) the block of N/4 frequency domain complex numbers includes a k-th frequency domain complex 
number, and 

(3) the pre-processing means includes: 



18 



EP 0 535 893 A2 



(a) a first deriving means for deriving, form the block of N/2 spectral coefficients, a block of first 
intermediate numbers, the first intermediate numbers including a /c-th, a 2/c-th, and a (2k + 1)-th 
first intermediate number, the first deriving means deriving the first intermediate numbers by: 

(i) selecting the 2/c-th spectral coefficient as the /c-th first intermediate number for each integral 
value of k between 0 and N/4 - 1, and 

(ii) reversing the sign of the (N - 1 - 2/c)-th spectral coefficient to provide a sign-reversed (N - 1 
- 2/c)-th spectral coefficient, and providing the sign-reversed (N- 1 - 2/c)-th spectral coefficient 
as the /c-th first intermediate number for each integral value of k between N/4 to NI2 - 1 , and 

(b) a second deriving means for deriving, from the block of N/2 first intermediate numbers, the block 
of N/4 frequency domain complex numbers, the frequency domain complex numbers having a real 
part and an imaginary part, and including a /c-th frequency domain complex number, the second 
deriving means deriving the block of N/4 frequency domain complex numbers by: 

(i) selecting the 2/c-th first intermediate number as a real part of a /c-th complex number, 

(ii) selecting the (2/c + 1)-th first intermediate number as an imaginary part of the /c-th complex 
number, and 

(iii) multiplying the /c-th complex number by an exponential of (-4njklN) to provide the /c-th fre- 
quency domain complex number, 

for each integral value of k between 0 and N/4-1. 

The inverse transform processing apparatus of claim 5, wherein, 

(1) block of N/2 processed samples includes an n-th processed sample, where n in an integer between 
0 and A//2- Land 

(2) the post- processing means includes: 

(a) a means for expanding the block of N/4 time domain complex numbers based on periodicity to 
provide a block of 3N/4 complex numbers, the block of 3N/4 complex numbers including an n-th 
complex number and a (N/4 - 1 - n)-th complex number, 

(b) a means for deriving a (N/4 - 1 - n)-th complex conjugate from the (N/4- 1 -n)-th complex number, 

(c) a means for 

(i) summing the (N/4 - 1 - n)-th complex conjugate and the n-th complex number to provide a 
first sum, and 

(ii) dividing the first sum by two to provide an n-th third intermediate complex number, 

(d) a means for 

(i) subtracting the (N/4 - 1 - n)-th complex conjugate from the n-th complex number to provide 
a difference, 

(ii) dividing the difference by 2j to provide a quotient, where j is the square root of minus one, 
and 

(iii) multiplying the quotient by the exponential of (-2nj(2n + 1)/N) to provide an n-th fourth in- 
termediate complex number, and 

(e) a means for 

(i) summing the n-th third intermediate complex number and the n-th fourth intermediate com- 
plex number to provide an n-th fifth intermediate complex number, 

(ii) multiplying the n-th fifth intermediate complex number by an exponential of (-nj(2n + 1)/(2N)) 
to provide an n-th sixth intermediate complex number having a real part and an imaginary part, 
and 

(iii) providing the real part of the n-th sixth intermediate complex number as the n-th processed 
sample. 

The inverse transform processing apparatus of claim 5, wherein 

(1) the block of N/2 processed samples includes an (n + N/4)-th, an (3N/4 - 1 - n)-th, and an (n - 3N/4)- 
th processed sample, each processed sample having a sign, 

(2) the block of N samples in the time domain includes an n-th sample in the time domain, and 

(3) the means for multiplying the block of N/2 processed samples in the time domain by an inverse 
transform window includes: 

(a) a deriving means for deriving from the block of N/2 processed samples, a block of N intermediate 
processed samples, the block of N intermediate processed samples including an n-th intermediate 
processed sample, the driving means deriving the block of N intermediate processed samples by: 

(i) selecting the (n + N/4)-th processed sample as the n-th intermediate processed sample for 

each integral value of n from 0 to N/4 - 1 , 
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(ii) reversing the sign of the (3N/4 - 1 - n)-th processed sample to provide a sign-reversed (3N/4 
• 1 - n)-th processed sample, and providing the sign-reversed (3N/4 - 1 - n)-th processed sample 
as the n-th intermediate processed sample for each integral value of n from N/4 to 3N/4 - 1, and 

(iii) reversing the sign of the (n - 3n)-th processed sample to provide a sign-reversed (n - 3N/4)- 
th processed sample, and providing the sign-reversed (n - 3A//4>-th processed sample as the 
n-th intermediate processed sample for each integral value of n from 3N/4 to N - 1; and 

(b) a means for multiplying the n-th intermediate processed sample by an inverse transform window 
for the n-th intermediate processed sample, for each integral value of n form 0 to N - 1, to provide 
the n-th sample in the time domain. 

A method for performing a forward modified discrete cosine transform on a block of N samples of a digital 
input signal in the time domain to provide a block of N/2 spectral coefficients in the frequency domain, 
the method comprising the steps of: 

(1) multiplying the block of N samples by a forward transform window to provide a block of N forward 
transform window processed samples; and 

(2) transforming the block of N forward transform window processed samples to provide a block of N/2 
spectral coefficients, the step of transforming the block of N forward transform window processed sam- 
ples including the steps of: 

(a) preprocessing the block of N forward transform window processed samples to provide a block 
of A//4 complex numbers, 

(b) performing processing equivalent to a Fast Fourier Transform on the block of N/4 complex num- 
bers to provide a transformed block of N/4 complex numbers; and 

(c) post-processing the transformed block of complex numbers to provide the block the N/2 spectral 
coefficients. 

The method for performing a forward modified discrete cosine transform of claim 9, wherein 

(1) the block of N forward transform window processed samples includes an n-th sample, a (n + 3N)- 
th sample, the (n + 3N/4>-th sample having a sign, and a (n - A//4)-th sample; and 

(2) the step of preprocessing the block of N forward transform window processed samples includes 
the steps of: 

(a) deriving, from the block of N forward transform window processed samples, a block of first in- 
termediate samples, the block of first intermediate samples including an n-th a (N- 1 - 2n)-th, and 
a 2n-th first intermediate sample, the block of first intermediate samples being derived by steps 
including: 

(i) reversing the sign of the (n + 3N/4)-th sample to provide a sign-reversed (n + 3N/4)-th sample, 
and providing the sign-reversed (n + 3N/4)-th sample as the n-th first intermediate sample, for 
each integral value of n from 0 to N/4 - 1 , and 

(ii) providing the (n - A//4)-th sample as the nth first intermediate sample, for each integral value 
of n from N/4 to N - 1 , 

(b) deriving, from the block of first intermediate samples, a block of second intermediate samples, 
the block of second intermediate samples including an n-th, a 2n-th, and a (2n + 1)-th second in- 
termediate sample, the block of second intermediate samples being derived by steps including sub- 
tracting, for each integral value of n from 0 to N/2 - 1, the (N - 1 - 2n)-th first intermediate sample 
from the 2n-th first intermediate sample to provide the n-th second intermediate sample, 

(c) deriving, from the block of second intermediate samples, the block of N/4 complex numbers, the 
complex numbers in the block of N/4 complex numbers each having a real part and an imaginary 
part, and including an n-th complex number, the block of N/4 complex numbers being derived by 
steps including: 

(i) selecting the2n-th second intermediate sample as a real part of an n-th intermediate complex 
number, 

(ii) selecting the (2N + 1)-th second intermediate sample as an imaginary part of the n-th inter- 
mediate complex number, and 

(iii) multiplying the n-th intermediate complex number by the exponent of (-4rc/n/N) to provide 
the n-th complex number in the block of N/4 complex numbers, 

for each integral value of n from 0 to N/4 - 1 . 

The method for performing a forward modified discrete cosine transform according to claim 9, wherein 
(1) the block of N/2 spectral coefficients includes a fc-th spectral coefficient, where k is an integer be- 



20 



EP 0 535 893 A2 



tween 0 and N/2 - 1, and 

(2) the step of post-processing the transformed block of complex numbers includes the steps of: 

(a) expanding the transformed block of N/4 complex numbers based on periodicity to provide a block 
of 3A//4 complex numbers, the block of 3A//4 complex numbers including a k-th complex number 

5 and a (N/4 - 1 - kyth complex number, 

(b) deriving a (N/4 - 1 - k)-th complex conjugate from the (NI4 - 1 - k)-th complex number, 

(c) deriving a k-th third intermediate complex number by 

(i) summing the (N/4 - 1 - k)-th complex conjugate and the k-th complex number to provide a 
sum, and 

w (ii) dividing the sum by two to provide the /c-th third intermediate complex number, 

(d) deriving a k-th fourth intermediate complex number by 

(i) subtracting the (N/4 - 1 - k)-th complex conjugate from the k-th complex number to provide 
a difference, 

(ii) dividing the difference by 2j to provide a quotient, where j is the square root of minus one, 
15 and 

(iii) multiplying the quotient by the exponential of (-2itj(2k + 1)/A/) to provide the /c-th fourth in- 
termediate complex number, 

(e) deriving the k-th spectral coefficient by 

(i) summing the k-th third intermediate complex number and the k-th fourth intermediate com- 
20 plex number to provide a k-th fifth intermediate complex number, 

(ii) multiplying the k-th fifth intermediate complex number by an exponential of (-nj(2k + \)/(2N) 
to provide a k-th sixth intermediate complex number having a real part and an imaginary part, 
and 

(iii) providing the real part of the k-th sixth intermediate complex number sum as the k-th spec- 
25 tral coefficient. 

12. The method of performing a forward modified discrete cosine transform according to claim 9, wherein, 

(1) the block of N/2 spectral coefficients includes a k-th spectral coefficient, where k is an integer be- 
tween 0 and N/2 - 1 , and 

30 (2) step of post-processing the transformed block of complex numbers includes the steps of: 

(a) expanding the transformed block of N/4 complex numbers based on periodicity to provide a block 
of 3A//4 complex numbers, the complex numbers in the block of 3A//4 complex numbers having a 
real part and an imaginary part and including a k-th complex number and a (N/4 - 1 - k)-th complex 
number; 

35 (b) determining: 

(i) a sum, the sum being one half of a sum of the cosine of (2n(2k + 1)/4W) and sine of (10n(2/r 
+1)/4N), 

(ii) a first difference, the first difference being one half of a difference between the cosine of 
(2n(2k + 1)/4/v) and the sine of (10n(2/c + 1)/4W), and 

40 (iii) a second difference, the second difference being one half of a difference between and sine 

of (10n(2/c + 1)/4W) and the cosine of (2n(2k + 1)/4N); 
(c) multiplying: 

(i) the real part of the k-th complex number by the first difference to provide a first intermediate 
number, 

45 (ii) the real part of the (N/4 - 1 - k)-th complex number by the sum, to provide a second inter- 

mediate number, 

(iii) the imaginary part of the k-th complex number by the sum to provide a third intermediate 
number, and 

(iv) the imaginary part of the (N/4 - 1 - k)-th complex number by the second difference to provide 
so a fourth intermediate number; and 

(d) summing the first, the second, the third, and the fourth intermediate numbers to provide the 
k-th spectral coefficient. 

13. A method of performing an inverse modified discrete transform on a block of N/2 spectral coefficients in 
55 the frequency domain to provide a block of N samples in the time domain as a digital output signal, the 

method comprising the steps of: 

(1) inversely transforming the block of N/2 spectral coefficients to provide a block of N/2 processed 
samples in the time domain, the step of inversely transforming the block of N/2 spectral coefficients 
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including the steps of: 

(a) pre-processing the spectral coefficients in the block of N/2 spectral coefficients to provide a 
block of N/4 frequency domain complex numbers, 

(b) performing processing equivalent to a Fast Fourier Transform on the block of N/4 frequency do- 
main complex numbers to provide a block of N/4 time domain complex numbers, and 

(c) post-processing the block of N/4 time domain complex numbers to provide the block of N/2 proc- 
essed samples in the time domain; and 

(2) re-arraying and multiplying the block of N/2 processed samples in the time domain by an inverse 
transform window to provide the block of N samples in the time domain. 

14. The method for performing an inverse modified discrete transform processing of claim 13, wherein, 

(1) the block of N/2 spectral coefficients includes a 2k- th, and a (N - 1 - 2/r)-th spectral coefficient, 
each having a sign, 

(2) the block of N/4 frequency domain complex numbers includes a k-th frequency domain complex 
number, and 

(3) the step of pre-processing the spectral coefficients in the block of N/2 spectral coefficients to pro- 
vide a block of N/4 frequency domain complex numbers includes the steps of: 

(a) deriving, from the block of N/2 spectral coefficients, a block of first intermediate numbers, the 
first intermediate numbers including a fc-th, a 2/r-th, and a (2k* 1)-th first intermediate number, the 
first intermediate numbers being derived by steps including: 

(i) selecting the 2k-th spectral coefficient as the fc-th first intermediate number for each integral 
value of k between 0 and N/4 - 1 , and 

(ii) reversing the sign of the (N- 1 -2/c)-th spectral coefficient to provide a sign-reversed (N - 1 
- 2/f)-th spectral coefficient, and providing the sign-reversed (N - 1 - 2/c)-th spectral coefficient 
as the k-th first intermediate number for each integral value of k between N/4 to N/2 - 1 , 

(b) deriving, from the block of N/2 first intermediate numbers, the block of N/4 frequency domain 
complex numbers, the frequency domain complex numbers having a real part and an imaginary 
part, and including a k-th frequency domain complex number, the block the frequency domain com- 
plex numbers being derived by steps including: 

(i) selecting the 2k-th first intermediate number as a real part of a k-th complex number, and 

(ii) selecting the (-2k + 1 )-th first intermediate number as an imaginary part of the k-th frequency 
domain complex number, and 

(iii) multiplying the k-th complex number by an exponential of (-4njk/N) to provide the k-th fre- 
quency domain complex number, 

for each integral value of k between 0 and N/4-1. 

15. The method of performing an inverse modified discrete cosine transform according to claim 13, wherein, 

(1) block of N/2 processed samples includes an n-th processed sample, where n is an integer between 
0 and N- 1, and 

(2) the step of pre-processing the block of N/4 time domain complex numbers includes the steps of: 

(a) expanding the block of N/4 time domain complex numbers based on periodicity to provide a block 
of 3A//4 complex numbers, the block of 3N/4 complex numbers including an n-th complex number 
and a (N/4 - 1 - n)-th complex number, 

(b) deriving a (N/4 - 1 - n)-th complex conjugate from the (N/4 - 1 - n)-th complex numbers, 

(c) summing the (N/4 - 1 - n)-th complex conjugate and the n-th complex number to provide a first 
sum, and dividing the first sum by two to provide an n-th third intermediate complex number, 

(d) deriving an n-th fourth intermediate complex number by 

(i) subtracting the (N/4 - 1 - n)-th complex conjugate from the n-th complex number to provide 
a difference, 

(ii) dividing the difference by 2j to provide a quotient, where j is the square root of minus one, 
and 

(iii) multiplying the quotient by an exponential of (-2nj(2n + 1)/N) to provide the n-th fourth in- 
termediate complex number, and 

(e) deriving a processed sample by 

(i) summing the n-th third intermediate complex number and the n-th fourth intermediate com- 
plex number to provide an n-th fifth intermediate complex number, 

(ii) multiplying the n-th fifth intermediate complex number by an exponential of (nj(2n + 1)/(2N)) 
to provide an n-th sixth intermediate complex number having a real part and an imaginary part, 
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and 

(iii) providing the real part of the n-th sixth intermediate complex number as the n-th processed 
sample. 

The method of performing an inverse modified discrete cosine transform according to claim 13, wherein 

(1) the block of N/2 processed samples includes an (n + N/4)-th, an (3N/4 - 1 - n>-th. and an (n - 3N/4)- 
th processed sample, each processed sample having a sign, 

(2) the block of N samples in the time domain includes an n-th sample in the time domain, and 

(3) the step of multiplying the block of N/2 processed samples in the time domain by an inverse trans- 
form window includes the steps of: 

(a) deriving, from the block of N/2 processed samples, a block of N intermediate processed samples, 
the block of N intermediate processed samples including an /7-th intermediate processed sample, 
the block of N intermediate processed samples being derived by steps including: 

(i) selecting the (n + A//4)-th processed sample selecting the (n + N/4)-th processed sample as 
the n-th intermediate processed sample for each integral value of n from 0 to A//4 - 1, 

(ii) reversing the sign of the (3A//4 - 1 - n)-th processed sample to provide a sign-reversed (3N/4 
- 1 - n)-th processed sample, and providing the sign-reversed (3A//4 - 1 - n)-th processed sample 
as the n-th intermediate processed sample for each integral value of n from N/4 to 3N/4 1 , and 

(iii) reversing the sign of the (n - 3N/4)-th processed sample to provide a sign-reversed (n - 3N/4)- 
th processed sample, and providing the (n - 3N/4)-th processed sample as the n-th intermediate 
processed sample for each integral value of n from 3 N/4 to N - 1 , and 

(b) multiplying the n-th intermediate processed sample by an inverse transform window for the n- 
th intermediate processed sample, for each integral value of n from 0 to N - 1, to provide the n-th 
sample in the time domain. 

A medium for storing a compressed digital signal comprising a block of N/2 spectral coefficients derived 
from a block of N samples of a digital input signal in the time domain by a modified discrete cosine trans- 
form performed according to a method comprising the steps of: 

(1) multiplying the block of N samples by a forward transform window to provide a block of N forward 
transform window processed samples; and 

(2) transforming the block of N forward transform window processed samples to provide a block of N/2 
spectral coefficients, the step of transforming the block of N forward transform window processed sam- 
ples including the steps of: 

(a) preprocessing the block of N forward transform window processed samples to provide a block 
of N/4 complex numbers, 

(b) performing processing equivalent to a Fast Fourier Transform on the block of N/4 complex num- 
bers to provide a transformed block of N/4 complex numbers; and 

(c) post-processing the transformed block of complex numbers to provide the block of N/2 spectral 
coefficients. 

18. The medium of claim 17, wherein, in the method of performing the forward modified discrete cosine trans- 
form, 

(1) the block of N forward transform window processed samples includes an n-th sample, a (n + 3N/4)- 
th sample, the (n + 3A//4)-th sample having a sign, and a (n - N/4)-th sample; and 

(2) the step of preprocessing the block of N forward transform window processed samples includes 
the steps of: 

(a) deriving, from the block of N forward transform window processed samples, a block of first in- 
termediate samples, the block of first intermediate samples including an n-th a (A/ - 1 - 2n)-th, and 
a 2n-th first intermediate sample, the block of first intermediate samples being derived by steps 
including: 

(i) reversing the sign of the (n + 3N/4)-th sample to provide a sign-reversed (n + 3N/4)-th sample, 
and providing the sign-reversed (n + 3N/4)-th sample as the n-th first intermediate sample, for 
each integral value of n from 0 to N/4 - 1, and 

(ii) providing the (n - A//4)-th sample as the nth first intermediate sample, for each integral value 
of n from N/4 to N - 1 , 

(b) deriving, from the block of first intermediate samples, a block of second intermediate samples, 
the block of second intermediate samples including an n-th, a 2n-th, and a (2n + 1)-th second in- 
termediate sample, the block of second intermediate samples being derived by steps including sub- 
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trading, for each integral value of n from 0 to N/2 - 1, the (N - 1 - 2n)-th first intermediate sample 
from the 2/7-th first intermediate sample to provide the n-th second intermediate sample, 
(c) deriving, from the block of second intermediate samples, the block of N/4 complex numbers, the 
complex numbers in the block of N/4 complex numbers each having a real part and an imaginary 
part and including an n-th complex number, the block of A//4 complex numbers being derived by 
steps including: 

(i) selecting the 2n-th second intermediate sample as a real part of an n-th intermediate complex num- 
ber, 

(ii) selecting the (2n + 1)-th second intermediate sample as an imaginary part of the n-th intermediate 
complex number, and 

(iii) multiplying the n-th intermediate complex number by the exponent of (-4rt/n/N) to provide the n-th 
complex number in the block of N/4 complex numbers, 

for each integral value of n from 0 top N/4-1 . 

19. The medium of claim 17, wherein, in the method of performing the forward modified discrete cosine trans- 
form, 

(1) the block of N/2 spectral coefficients includes a /c-th spectral coefficient, where k is an integer be- 
tween 0 and N/2 - 1, and 

(2) the step of post-processing the transformed block of complex numbers includes the steps of: 

(a) expanding the transformed block of N/4 complex numbers based on periodicity to provide a block 
of N/2 complex numbers, the block of N/2 complex numbers including a Ac-th complex number and 
a (N/4 - 1 -/c)-th complex number, 

(b) deriving a (N/4 - 1 - /c)-th complex conjugate from the (N/4 - 1 - /c)-th complex number, 

(c) deriving a /c-th third intermediate complex number by 

(i) summing the (N/4 - 1 - k)-th complex conjugate and the /c-th complex number to provide a 
sum, and 

(ii) dividing the sum by two to provide the Ac-th third intermediate complex number, 

(d) deriving a /c-th fourth intermediate complex number by 

(i) subtracting the (N/4 - 1 - /c)-th complex conjugate from the /c-th complex number to provide 
a difference, 

(ii) dividing the difference by 2j to provide a quotient, where j is the square root of minus one, 
and 

(iii) multiplying the quotient by the exponential of (-2nj(2k + 1)/N) to provide the /c-th fourth in- 
termediate complex number, 

(e) deriving the /c-th spectral coefficient by 

(i) summing the /c-th third intermediate complex number and the /c-th fourth intermediate com- 
plex number to provide a /c-th fifth intermediate complex number, 

(ii) multiplying the /c-th fifth intermediate complex number by an exponential of (-Ji/(2k + 1)/(2N)) 
to provide a /c-th sixth intermediate complex number having a real part and an imaginary part, 
and 

(iii) providing the real part of the /c-th sixth intermediate complex number sum as the /c-th spec- 
tral coefficient. 

20. The medium of claim 17, wherein, in the method of performing the forward modified discrete cosine trans- 
form, 

(1) the block of N/2 spectral coefficients includes a /c-th spectral coefficient, where k is an integer be- 
tween 0 and N/2-1, and 

(2) step of post-processing the transformed block of complex numbers includes the steps of: 

(a) expanding the transformed block of N/4 complex numbers based on periodicity to provide a block 
of 3N/4 complex numbers, the complex numbers in the block of 3N/4 complex numbers having a 
real part and an imaginary part and including a /c-th complex number and a (N/4 - 1 - /c)-th complex 
number; 

(b) determining: 

(i) a sum, the sum being one half of a sum of the cosine of (2jt(2/c + 1)/4N) and sine of (10ir(2A 
+ 1)/4N), 

(ii) a first difference, the first difference being one half of a difference between the cosine of 
(2n(2/c + 1)/4N) and the sine of (10n(2/c + 1)/4N), and 

(iii) a second difference, the second difference being one half of a difference between and sine 
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of (10it(2/c + 1)/4W) and the cosine of {2n(2k + 1)/4N); 

(c) multiplying: 

(i) the real part of the fc-th complex number by the first difference to provide a first intermediate 
number, 

5 (ii) the real part of the (N/4 - 1 - /c)-th complex number by the sum, to provide a second inter- 

mediate number, 

(iii) the imaginary part of the fr-th complex number by the sum to provide a third intermediate 
number, and 

(iv) the imaginary part of the (A//4 - 1 - /c)-th complex number by the second difference to provide 
10 a fourth intermediate number; and 

(d) summing the first, the second, the third, and the fourth intermediate numbers to provide the k- 
th spectral coefficient. 

21. The medium of claim 17, wherein the medium is an optical disc. 

15 

22. The medium of claim 17, wherein the medium comprises a semiconductor memory. 

23. The medium of claim 17, wherein the medium comprises a magnetic tape. 
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